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Abstract 

Quantum and thermal fluctuations of an irrotational fluid are studied across the transi- 
tion regime connecting a protoinfiationary phase of decelerated expansion to an accelerated 
epoch driven by a single inflaton field. The protoinfiationary inhomogeneities are suppressed 
when the transition to the slow roll phase occurs sharply over space-like hypersurfaces of 
constant energy density. If the transition is delayed, the interaction of the quasi-normal 
modes related, asymptotically, to fluid phonons and inflaton quanta leads to an enhance- 
ment of curvature perturbations. It is shown that the dynamics of the fluctuations across 
the protoinfiationary boundaries is determined by the monotonicity properties of the pump 
fields controlling the energy transfer between the background geometry and the quasi-normal 
modes of the fluctuations. After corroborating the analytical arguments with explicit nu- 
merical examples, general lessons are drawn on the classification of the protoinfiationary 
transition. 
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1 Introduction 



In the conventional lore, the large-scale temperature and polarization anisotropies of the 
Cosmic Microwave Background are caused by curvature inhomogeneities with typical wave- 
lengths exceeding the Hubble radius at the time of matter radiation equality [I1I2]- A nearly 
fiat spectrum of Gaussian fluctuations of the spatial curvature naturally arises from the 
quantum inhomogeneities of a single infiaton field evolving during a quasi-de Sitter stage of 
expansion. Although the simplest scenario is consistent with the observational signatures, 
different sets of initial conditions have been explored through the years. 

Initial states different from the vacuum can modify the temperature and polarization 
anisotropies at large scales. This general approach has been scrutinized along various per- 
spectives (see, e.g. [21 HI O [6], [71 El O [101 E])- Temperature-dependent phase transitions 
[21 [1] lead to an initial thermal state for the metric perturbations [3 [6l [71 [9l [10] • If the initial 
state is not thermal (but it is not the vacuum either), curvature phonons can be similarly 
produced via stimulated emission. Second-order correlation effects of the scalar and tensor 
fluctuations of the geometry can be used to explore the statistical properties of the initial 
quantum state [8] by applying the tenets of Hanbury-Brown-Twiss interferometry [12] which 
is employed, in quantum optics, to infer the bunching properties of visible light. 

The modifications of the initial state are subjected to a number of constraints all orig- 
inating, directly or indirectly, from the comparison between the energetic content of the 
initial fluctuations and the energy density of the background geometry. The criterion for the 
avoidance of severe backreaction effects is not unique. Single field quasi-de Sitter inflationary 
models with general initial states of primordial quantum fluctuations have been examined in 
[HI [H] with the aim of deriving constraints from the study of higher-order correlation func- 
tions and from the renormalizability of the energy-momentum tensor of the fluctuations. 
It is equally plausible to demand that the energy-momentum pseudo-tensor of the scalar 
and tensor fluctuations does not exceed the energy density and pressure of the background 
geometry, as argued in [T3] . 

In the present paper a complementary and novel approach to the problem of the initial 
conditions of cosmological perturbations is pursued. The ever expanding inflationary back- 
grounds are geodesically incomplete and inflation cannot be eternal in the past. Thus it is 
legitimate to suppose the existence of a protoinflationary phase where the dynamics of the 
background was not yet accelerated. As the terminology indicates, the purpose here is not 
to test the universality of inflation given a set of arbitrary and widely different choices of 
the preinflationary dynamics. While under certain conditions inflation can be dynamically 
realized, it is not eternal in the past either. The modest purpose here is not to challenge 
inflation but to analyze the initial conditions of large-scale inhomogeneities in an improved 
dynamical framework. Just to avoid potential misunderstandings, it should be clear that the 
transition from a deceleration to acceleration has nothing to do with the so-called bouncing 
behaviour where the background passes from contraction to expansion or vice versa (see, e.g. 
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[2] and references therein). In the present framework the universe will always be expanding 
even during the protoinflationary phase. 

The approach suggested here is pragmatic and the attention is focused on single field in- 
flationary models leading to a nearly fiat spectrum of curvature inhomogeneities [15]. During 
the protoinflationary phase of decelerated dynamics the energy momentum tensor is domi- 
nated by a single perfect fluid. The analysis can be generalized to include various inflaton 
fields and protoinflationary fluids but this will not be the primary goal of this investigation. 
Unlike the standard scenario, during the protoinflationary phase the seeds of curvature inho- 
mogeneities are fluid phonons, i.e. the quantum excitations of an irrotational and relativistic 
fluid discussed by Lukash [TB] (see also [T71 HB]) right after one of the first formulations of 
inflationary dynamics [T^]. The whole irrotational system can be reduced to a single (de- 
coupled) normal mode which is promoted to field operator in case the initial fluctuations 
are required to minimize the quantum Hamiltonian of the phonons [16]. The canonical nor- 
mal mode identified in [16] is invariant under infinitesimal coordinate transformations as 
required in the context of the Bardeen formalism [20] (see also [IT]). The subsequent analy- 
ses of Refs. [21] and [22] follow the same logic of [16] but in the case of scalar field matter; the 
normal modes identified in [161 ISl |22] coincide with the (rescaled) curvature perturbations 
on comoving orthogonal hypersurfaces [231 [23] (see also the beginning of section [3|). 

The fluid phonons can be treated quantum mechanically but the initial state does not 
need to be the vacuum: if the protoinflationary phase is dominated by radiation, the fluid 
phonons are more likely to follow a Bose-Einstein distribution as it happens in inflationary 
models based on temperature-dependent phase transitions [H [5l [6]. When the protoinfla- 
tionary inhomogeneities are suppressed across the boundary, the initial normalization of 
curvature perturbations is set, most likely, by the quantum mechanical fluctuations gen- 
erated during the inflationary phase and, depending on the duration of inflation, by the 
stimulatated emission from the protoinflationary relics. The conditions for the suppression 
or for the enhancement of protoinflationary curvature perturbations are related to the evolu- 
tion of the pump fields controlling the transfer between the energy density of the background 
and the quasi-normal modes of the system. 

The present paper is organized as follows. After introducing the governing equations, in 
section |21 the quasi-normal mode of the system are derived in section [31 The fate of the 
large-scale curvature perturbations across the protoinflationary transition is investigated in 
section [H where the normalization of the fluctuations during the protoinflationary phase is 
also discussed. In section \5\ the nature of the transition is clarified in terms of the mono- 
tonicity properties of the pump fields accounting for the contribution of the inflaton and 
of the protoinflationary fluid to the total curvature perturbations. Explicit analytical and 
numerical examples are used to draw some general lessons on the dynamical features of the 
protoinflationary transition. Section ^ contains some concluding remarks. The derivation of 
the coupled evolution equations of the quasi-normal modes related, asymptotically, to fluid 
phonons and inflaton quanta is reported in the appendix. 
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2 Governing equations 



2.1 General consideration 



The minimal set of assumptions characterizing the framework of the present investigation 
stipulates that the four-dimensional geometry is determined by the Einstein equations, sup- 
plemented by the conservation equations accounting for the dynamics of the inflaton and of 
the protoinflationary source^: 



dV_ 

dip 



0, 



0, 



g°'^UaUi3 = 1, 



(2.1) 

(2.2) 
(2.3) 



where T^{y^) and T^{ppr, Ppr) are, respectively, the energy-momentum tensors of the inflaton 
field (fi and of the protoinflationary fluid: 



"TaiPpr, Ppr) 



(Ppr ~l~ Ppr) UaU^ Ppr S^. 



(2.4) 
(2.5) 



The subscripts in the energy density and pressure remind of the protoinflationary origin of 
the fluid variables. In a conformally flat background metric of the type ^^.^ = a^{T)r]ai3 
(where a(r) is the scale factor in conformal time and rja/s is the Minkowski metric), Eqs. 
(12. ip . (12. 2 p and (12. 3p lead to a set of four equations 



8nG 



2 V 2 

a Ppr + — + V a 



dV 

^" + 2n^' + —c? = 0, 

p' + Snippr+Ppr) = 0, 



(2.6) 

(2.7) 
(2.8) 
(2.9) 



which are not all independent and whose specific form is dictated by the fluid content of the 
primordial plasma. In Eqs. fl2.6l) - (l2.9p the prime denotes a derivation with respect to the 
conformal time coordinate r; furthermore H = (In a)'. The connection between H and the 
Hubble parameter is H = Ti/a. The effective energy and pressure densities of ip are given 
by 



./2 



Pv 



2a 



,+V{p>), 



P 



/2 

2a? 



(2.10) 



^ Greek indices run from to 3; the signature of the metric gai3 is mostly minus and Va denote the 
covariant derivative with respect to gap ■ 
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Note that Eq. (12. 8p is equivalent to 



p;.3«(p,.p.)^o, 4^J^i.|^(|)« (2.11) 

2.2 Uniform curvature gauge 

The most general scalar fluctuation of the four-dimensional metric is parametrized by four 
different functions whose number can be eventually reduced by specifying (either completely 
or partially) the coordinate system: 

45-00 = 2a^0, 5^gij = 2a^{ip5ij - didjE), S^goi = -o^diB, (2.12) 

where Sg denotes the scalar mode of the corresponding tensor component; the full metric 
(i.e. background plus inhomogeneities) is given, in these notations, by ^^/^(a;, r) = ^^^(r) + 
Ssgai3{x,T) where, as already mentioned prior to Eqs. fl2.6l) - (l2.9p Qapi'T) = o?'{T)rjaj3- For 
infinitesimal coordinate shifts r — )■ r = r + eo and — )■ = + 9*e the functions 0(x, r), 
B{x,t), iP{x,t) and E{x,t) introduced in Eq. (I2.12p transform as: 

0^0 = 0-7^eo-e[), ^^iJ = ^ + 'Heo, (2.13) 
B = B + eo- e', E = E - e. (2.14) 

In the uniform curvature gauge two out of the four functions of Eq. ( 12. 12^ are set to zero 
[25]: 

E = 0, ip = 0. (2.15) 

Starting from a gauge where E and ip do not vanish, the perturbed line element can always be 
brought in the form fl2.15p by demanding E = and ip = in Eqs. f l2.13p and fl2.14p . More 
specifically, if E ^ and ip ^ 0, the uniform curvature gauge condition can be recovered by 
fixing the gauge parameters as e = E and eo = —ip/T-L. This choice guarantees that, in the 
transformed coordinate system, ip = E = 0. 

The gauge condition of Eq. fl2.15p implies that the fluctuations of the spatial curvature 
vanish but, in this case, the perturbed metric also contains off-diagonal elements. With the 
condition fl2.15p the gauge freedom is totally fixed without the need of further conditions: 
because of this property the functions 0(x, r) and B{x, r) bear an extremely simple relation 
to one of the conventional sets of gauge-invariant variables, as it will be shown later in this 
section. Finally, the off-diagonal coordinate system will prove very useful in section [3] and 
in appendix |A] for the analysis of the coupled system of quasi-normal modes. In the gauge 
fl2.15p the inhomogeneities of the energy-momentum tensors T^{ip) and XiiPpr, Ppr ) ^1*6 f 
respectively, 

<5sT° = ^(-<Pip" + + xV), = --,v'd\ - ^d'B, (2.16) 

\ dip J 
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dV 



6p. 



pri 



5sV = -SPpr Si, 5J'q = {ppr + ppr)v\ 



(2.17) 

(2.18) 



where x{^: and f * denote, respectively, the fluctuations of the inflaton ip and the three- 
velocity field in the gauge (12.151) . The fluctuations of the Einstein tensor = R'^ — R^^^/'^ 
in the gauge (12.151) are instead: 

2 r .._o„ .... ,1 ^2.19) 



B' + 2nB] 5^ 



—didUB' + 2nB + 



-H(f)+{H' -H^)B 



(2.20) 
(2.21) 



The combination of Eqs. (l2A6|) -( l2Tri) and (12J8|) with Eqs. f l239|) . (I2:20|) and dMI]) implies 
that the (00) and (Oi) components of the perturbed Einstein equations with mixed indices 
become: 



r) = div\ The variables 6p^ and 6py, c< 
energy density and of the pressure of the inflaton field: 



'2v72; 



(2.22) 

(2.23) 



where 6pr{x,T) = div\ The variables (5p<^ and 6p^ correspond to the fluctuations of the 



Sp. 



f 9 



(2.24) 
(2.25) 



To avoid lengthy notations we wrote 6p^ (instead of SsP^^), Sppr (instead of SsPpr) and similarly 
for the corresponding pressures; this notation is fully justified and unambiguous once the 
scalar nature of the fluctuations has been established, as specified by the general formulae 
written above. Bearing in mind this specification, the (ij) component of the perturbed 
Einstein equations reads: 

1. 



{n^ + 2n')(f) - Hcj)' - -V\(f) + B' + 2HB) 



SI 



+-did^ 
2 



B' + 2nB 



AnGa' 



5ppr + 6p^ 



6Um 



(2.26) 



The separation of the traceless part from the trace in Eq. (I2.26P leads to two independent 
relations: 

{n^ + 2'H')(j) + W + ^V'(0 + B' + 2'HB) = 4TTGa\6ppr + 6p^), (2.27) 
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di&>[(j) + B' + 2UB] - + B' + 2nB]5{ = SvrGa^n]. 

3 



(2.28) 



If the anisotropic stresses 11^ is neglected, Eq. f l2.26p is the equivalent to the following pair 
of conditions: 

{H^ + 2H')(f) + H(f)' = AwGa^Sppr + Sp^), (2.29) 
+ 5' + 27/5 = 0. (2.30) 

The evolution equation for the perturbed inflaton is: 

d'^V dV 

X" + 2'Hx' - V\ + Q^ci\ + 2(f)— - <^V' - ^p'V^B = 0. (2.31) 

Finally, by perturbing the covariant conservation equation of the fluid energy-momentum 
tensor the evolution equation for the density fluctuation is: 

Sppr + iP + P)0pr + 3H(5pp, + 5ppr) = 0, (2.32) 

while the equation for the three-divergence of the velocity becomes: 

(V + V^i?)' + [Pp'- + '^^Ppr + Ppr)] (g^^ ^ ^2^) ^ V%p. ^ ^ ^ ^2.33) 



Equations (12.311) and (12.331) represent the staring point for the derivation of the quasi-normal 
modes of the system. 



2.3 Gauge-invariant observables 

The coordinate system defined by Eq. (I2.15P completely fixes the gauge freedom without 
the need of further subsidiary conditions. The absence of spurious gauge modes is then 
guaranteed as it happens for other choices of coordinates removing completely the gauge 
freedom such as the conformally Newtonian gauge. As a consequence, the perturbation 
variables defined in the gauge (I2.15P bear a simple relation to the various gauge-invariant 
observables. The specific connection between the degrees of freedom defined in the uniform 
curvature gauge and other common gauge-invariant combinations will now be outlined. 

The curvature perturbation on comoving orthogonal hypersurfaces (conventionally de- 
noted by TZ) and the Bardeen potential (conventionally denoted by \E') coincide, up to time- 
dependent functions, with 0(x, r) and B{x,t): 

^ = -nB, n = -j^^<l^. (2.34) 

The result of Eq. (12.340 can be easily derived from the explicit gauge transformation re- 
lating the uniform curvature hypersurfaces with the comoving orthogonal hyepersurfaces. 
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Conversely, from the customary expression of TZ in terms of the gauge-invariant Bardeen po- 
tentials $ and \1/ the result of Eq. fl2.34p can be cross-checked. The curvature perturbations 
on comoving orthogonal hypersurfaces are given by: 

n.-,-lM^. (2.35) 

But in terms of the variables introduced in Eq. fl2.12p the expression of $ and \1/ is: 

^ = (j)+(^B-E'y + n{B-E'), ^ = iIj -H{B ~ E'). (2.36) 

According to Eq. (IT^ . in the gauge f EJf^ $ = + 5' + and ^ = -liB. By 
inserting into Eq. (I2.35P the expressions for $ and \1/ written in the gauge fl2.15p the results 
of Eq. (12.341) are immediately recovered. The total density contrast on uniform curvature 
hypersurfaces can be expressed as 

C = -^H, 5pt = {6ppr + 6p^), pt = pip + ppr- (2.37) 
Pt 

From Eq. (I2.22p recalling Eqs. (I2.34p and (I2.37P we can also obtain the relation between 7^, 
C and 

127rGa2(pt + Pt) 

It is relevant to remind that C, and TZ are often used interchangeably. This is justified provided 
the wavelengths under considerations are sufficiently larger than the Hubble radius at the 
corresponding time. Otherwise the two variables C, and TZ are physically different. 



3 Quasi- normal modes of the system 

The system of section [2] describing the evolution across the protoinflationary boundary has 
two asymptotic limits corresponding to the situation where one of the two components is 
either absent or dynamically negligible. If the protonflationary fluid and the infiaton are 
simultaneously present the evolution is characterized by a pair of (interacting) quasi-normal 
modes which are the generalization of the normal modes obtainable in the case of a single 
component. After swiftly summarizing what happens in the two asymptotic limits, the 
derivation of the quasi-normal modes of the whole system will be presented. The interested 
reader may also consult appendix |X] where some of the technical results involved in the 
derivation are collected. 

In the limit v''(r) — )■ and x(a;, r) — )■ the fluctuations of the infiaton energy density 
and of the infiaton pressure are both vanishing, i.e. dp^p = 6p^ = 0. Since 6ppr = c^Sppr, Eq. 
(I2.22P can be multiplied by and summed to Eq. (12.291) . After simple algebra the following 
result will be obtained: 

n(f)' + [H\l + 3cl) + 2n'](l) = -HclV^B. (3.1) 
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Introducing the variables a and 7 defined a; 



3. 



Eq. (13.ip becomes 



I A 2- (3-3) 

Differentiating botli sides of Eq. (13. 3 p witli respect to r, two kinds of terms (proportional 
to V^-B and to V^-B') will arise; using then Eq. (I2.29P to eliminate the terms proportional 
to V^B' and Eq. (13. 3p to eliminate the terms containing V^-B, a decoupled equation for 
is readily obtained. Recalling the simple relation between and TZ mentioned in (I2.34p the 
resulting equation becomes: 

z' 

7^;; + 2^7^;, - c^v'7^p, = o, (3.4) 

Zp^ 

where TZpr is the curvature perturbation on comoving orthogonal hypersurfaces and 



7 Zpr it Cs 

The function introduced in Eq. (13. 5p is actually the normal mode of the system obeying: 

- clV^q, - ^q, = 0. (3.6) 

Zpr 

By using the momentum constraint in the purely hydrodynamical case (i.e. ^p'{t) — )■ and 
x{x, r) — i- in Eq. (I2.23P ) the following chain of equalities holds 

A-nGa^^Ppr + ppr) 
^77 7 v. 

which also implies, always neglecting the inflaton, that V^T^pr- = T^i^pr + V^-B). The same 
steps leading to Eqs. (13. 4p and (13. 5 p can be applied to the case of scalar field matter in the 
absence of protoinflationary component (i.e. Ppr = Ppr = and 6ppr = Sppr = 9pr = 0). The 
evolution equation of the curvature perturbation will be given, in this case, by: 



v'Upr = -v'm = — '^';^^'^\ epr + v'B), (3.7) 



z', 



n'' + 2^n' - v'7^^ = 0, z^ = (3.8) 



Zip tt 



The separate use of the momentum constraint in the scalar field case (i.e. 6pr = in Eq. 
([223D) leads to the analog of Eq. ^l\f : 

= — = -f , (3.9) 

tt Zip 



•^The background-dependent functions and 7 mentioned in Eqs. p.ip - p.2p illustrate, for convenience, 
the notations Refs. [16j[l8]. In the rest of the paper it will be more practical to adopt a slightly different 
set of variables which are introduced in Eqs. p.4p and p.Sp . 
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where = ax = —z^ TZ^ obeys, from Eq. 



the following equation: 
0, 



(3.10) 



which is the analog of Eq. (13. 6 p holding in the absence of inflaton contribution. When dealing 
with the process of parametric amplification, the variables Zpr{t) and Zy,{t) are dubbed pump 
fields since they control the rate of energy transfer from the background to the fluctuations. 
This terminology is often used in quantum optics [12] (see also [S] ) and shall also be employed 
in the forthcoming considerations. 

The results obtained in Eqs. (I3.4I) - (I3.7I) assume the absence of the inflaton field. Con- 
versely Eqs. ( 13. Sp and ( 13. 9p assume the absence of the protoinflationary fluid. If the contribu- 
tions of the fluid and of the scalar field are simultaneously taken into account, the evolution 
equations of the resulting system can be reduced to a pair of coupled equations whose solu- 
tion gives directly the curvature fluctuations on comoving orthogonal hypersurfaces. Indeed, 
by keeping both contributions, Eq. (I2.23P implies 



4:71 G ip' 4:TTGa'^(ppr + Ppr) 

-X H —u, 



n 



n 



(3.11) 



where the notation {Op^ + B) = — V^m has been introduced. To simplify the discussion, 
we shall assume that the protoinflationary fluid is characterized by a constant barotropic 
index w so that Cs = y/w. The derivation of the coupled system of the quasi-normal modes 
is reported in the appendix |Al The final equations obeyed by qy and q^ are: 



- CsV^g^ + A,y{T)qy + By^{T)q^ + Cy^{T)q'^ = 0, 



(3.12) 
(3.13) 



The coefficients Ayy^r), By^ij) and C„^(r) depend on the conformal time coordinate and 
are given by the following expression^: 



Ayy{r) = - -^^ W- 7 - P^^^ L'\W - 1 



By^ir) 



H^{?>w + l) + 2H' 



2a\w + l)— + 2W{l-w)\ 

OLD ) 



^'^(w - 1) - 2(H^(3w + 1) + 2H') 



2'Hy/w 



w 



(3.14) 

(3.15) 
(3.16) 



"'in Eqs. p.l4p - p.l6p and p.l7p - (|3.19p natural Planckian units Mp = 1 are used. The same units are 
also used in the second part of appendix [A] where the explicit derivation of Eqs. p.l2p and (|3.13p is reported. 
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The coefficients A^^{t), B^y^r) and C^„(r) are instead: 



AHV dip n 



W 



1 



to 



yPpr + Ppr \ (3W - 1) (W - 1) I la^^^ 

9^9 



4?/ 



+ ^ 



w 



1 



2H '^\w 



(3.17) 

(3.18) 
(3.19) 



In the hmit Ppr — )■ and Ppr — )■ 0, the equation for coincides with the equation obeyed by 
z^plZ^ and reported in Eq. fl3.10p : this is not a surprise since, in the gauge f l2.15p and when 
Ppr = 0, = —z^ IZ^. In the opposite limit, i.e. ip' — )■ 0, the equation obeyed by g„ coincides 
with the equation obeyed by ZprlZpr (in the case Cg = w) and reported in Eq. (13.61) . In the 
general case, using the notations developed so far, the total curvature perturbations can be 
written as 

(3.20) 



From equation f l3.20p the curvature perturbations can also be computed and they are: 



n 



Zpr Cg 



Zpj.Cg + Z'^ 



-Qv - 



(3.21) 



It is immediate to show from Eqs. (13. 5p and (I3.2ip that TZ — t- TZp^ when 2;^ — )■ 0. In the same 
way when Zpr — )■ Eqs. (13. 9 p and (I3.2ip imply that 71 — )■ 7l^p. The quasi- normal modes qy 
and q^ describe, asymptotically, the excitations corresponding to fluid phonons and infiaton 
quanta. 



4 Across the protoinflationary transition 

The crudest model for the transition stipulates that prior to the onset of inflation the per- 
fect and irrotational fluid dominates the total energy density and pressure. In a slightly 
inhomogeneous space-time the energy-momentum tensor experiences a finite discontinuity 
on the hypersurface of constant energy density when the slow-roll dynamics starts off. This 
choice for the matching hypersurface is adopted in the case of post-inflationary transitions 
|26] and it is interesting to scrutinize its implications in describing the protoinflationary 
boundary. The logic of the sudden approximation is to assume that Zi^{t) — )■ during the 
protoinflationary phase. Conversely during the slow-roll epoch Zpr{t) — )• 0. The sudden 
approximation (together with the required continuity of the extrinsic curvature) implies the 
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suppression of the density contrast and of the metric fluctuation across the protoinflationary 
boundary. The potential limitations of the sudden approximation are scrutinized in section 
[5] where the protoinflationary dynamics is described in terms of a class of exact solutions of 
the background equations which will be presented later. 



4.1 The continuity of the extrinsic curvature 



If the stress tensor undergoes a finite discontinuity on a space-like hypersurface the inho- 
mogeneities are matched by requiring the continuity of the induced three metric and of the 
extrinsic curvature on that hypersurface. The extrinsic curvature is defined as 



1 

2iV 



V.iV, + V,iV, - 74 



(4.1) 



Recalling Eq. fl2.12p . in a generic coordinate system the lapse function, the shift vectors and 
three-metric jij are, respectively: 



;i + 2^)a^ 



-fij = a^(l - 2ij)Sij + 2a'^didjE. 



(4.2) 



Using Eq. ( 14. 2 p into Eq. (14. ip the covariant and mixed components of the extrinsic curvature 
read, to first order in the scalar metric perturbations. 



Kij{x, t) 



Kiix, t) 



—oHSij + a 



didj{B -E - 2HE) + (^A + 2H^/' + H^)5,j 



(4.3) 



The continuity of the background extrinsic curvature implies that across the transition the 
scale factor and the Hubble rate must be continuous. In cosmic time, a continuous form of 
the scale factor can be written when, for instance, the inflationary phase is characterized 
by a set of constant slow-roll parameters (see, e.g. Eq. (I4.37P of section H] for a general 
definition of the slow roll parameters). In this case we shall have that: 



Clpr (t) 



t Y 

a / t 



t < U 



+ 



(3 — a 



1/3 



t > U. 



(4.4) 



In Eq. (14. 4p the inflationary evolution is realized for /3 ^ 1. For some applications it is 
useful to recall the conformal time parametrization where the continuity of the scale factors 
across the protoinflationary boundary can be expressed as: 



apr(r) 
ai„/(r) 



P J Ti Ti 



a 



a r* 



-( = + l)r*<r< -r* 



-T* < T < — Ti, 



(4.5) 
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where a = a/{l — a) and (3 = (3/{(3 — 1). From Eq. fl4.4p the scale factor and its first 
derivative are continuous in t^,. By going in the conformal parametrization a(r) dr = dt the 
time coordinate r becomes negative and therefore the conformal time scale factor and its 
first derivative with respect to r are continuous in — r*. In the parametrization of Eq. (14.51) 
the infiationary regime occurs for /3 ^ 1 and /3 — )■ 1. 

From Eq. ( 14. 3 p the continuity of the inhomogeneous part of KI and Kij implies the 
separate continuity of the combinations 



0, 



'E] 



0, 



0, 



[B-E']± = 0, 



(4.6) 



where, following the general treatment of sudden transitions [26], the subscript ± denotes 
the jump of the corresponding quantity across the transition (i.e. [/]± = /+ — /-)• In the 
coordinate system where r is constant, the equation for the hypersurface of constant energy 
density becomes = 0. But since 



6pt 6p^ = 6pt - p[eo, 6p^ = 6ppr + (5p, 



if} 



(4.7) 



the condition Sp^ = implies eo = 6pt/p[. Recalling now the expressions for (p, ip, E and B 
stemming from Eqs. ( I2.13P and ( I2.14p . the conditions of Eq. ( 14. 6 p become 

= 0, \b-e'+^P' 

± 



pi 



Pt 



0, 



Pt 



0. 



(4i 



Equations (14.80 are general and can be studied in any gauge. According to Eq. (14. 8p . the con- 
tinuity of the extrinsic curvature in the gauge ( I2.15P implies that the following combinations 
must separately be continuous 



B + 



5pt 
Pi 



0, 



5pt 
Pi 



Pt 



± 

0. 



0, 



Since T-L is continuous, Eq. ( I4.10p reduces to 



Pi 



0. 



± 



The Hamiltonian constraint of Eq. (I2.22p can be written in the form: 

n 



Spt 



■[V'B + 3'H(j)], 



Equations (14. 9 p and (14. lip are then equivalent to the two conditions: 



0, 



0. 



(4.9) 
(4.10) 

(4.11) 

(4.12) 

(4.13) 
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But according to Eq. flXM]l we have that ^ = -HB and 7^ = -W^/iU'^ - 'H')(f)] thus 
the continuity of the scale factor and of "H imphes the continuity of TZ and \E' across the 
protoinflationary transition. The evolution will now be separately solved during the protoin- 
flationary phase and during the inflationary phase. The mathching conditions expressed by 
Eq. f l4.13p agree with former treatments [26j but in the coordinate system defined by Eq. 



4.2 Protoinflationary evolution 



During the protoinflationary phase and for — )■ the system reduces to the triplet of 
equations 



B + 2HB + 



\H J a 
2V'B H 



0, 



5, 



6ppr 
Ppr 



(4.14) 
(4.15) 
(4.16) 



where the overdot denotes a derivation with respect to the cosmic time coordinate and St 
is the total density contrast which is dominated, in this case, by the protoinflationary fluid. 
Equation ( I4.14p comes from Eq. f l2.30p : Eq. fl4.15p is Eq. (13. ip but written in the cosmic 
time coordinate; Eq. (I4.16P derives from Eq. ( I2.22p . 

To enforce a correct normalization on the solutions we proceed as follows. After promot- 
ing the normal mode q^{x,T) and its conjugate momentum to the status of field operators 
obeying canonical commutation relations at equal times, the mode expansion for g^,(x, r) 
becomes 



1 



(27r)3/2 

where [d^, at] = 5'^^^{k — p) and the mode function /^(r) obeys 

\/Vpr ~l~ Ppr 



(4.17) 



e ci - 



•^pr 



pr 



- 1/4 
V 



fk = 0, 



y 



^pr 



He, 



a 



?>{l-w) 
2(3w + 1)' 



(4.18) 
(4.19) 



In the case w = 1/3 (corresponding to a radiation fluid) z/ = 1/2; it is interesting to remark 
that, in the case w = (but with Cg 7^ 0), z/ = 3/2 allowing for a flat spectrum of phonons, 
as discussed in [16]. Equation (14. ISp can be solved exactly in terms of Hankel functions and 
the solution is: 



fk{T) 



-iTT{2u+l)/4 



(4.20) 
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where, as usual, we shall focus on the case Cg = y/w. If the initial conditions are not quantum 
mechanical but rather thermal, then the initial state will contain thermal phonons, i.e. 



(4.21) 



where T = aT denotes the comoving temperature while T is the physical temperature. The 
power spectrum of the fluid phonons can be easily determined from the two point function 
evaluated at equal times: 



(g„(f,r)g^(y,r)) 



sin k r 

kr ' 



\x-y\, 



2n' 



:\MT)fi2nk + l). 



(4.22) 



If kcs ^ T we have that — ?■ and the quantum mechanical initial conditions dominate; 
conversely if kc^ -C T the thermal initial conditions dominate against the quantum mechan- 
ical ones. Once the phonon spectrum is known the spectrum of curvature perturbations 
is ^ 

Vnik, t) = -A^|/,(r)|^ (271^ + 1). (4.23) 
27r^ z^.^ 

From the spectrum of curvature phonons it is elementary to derive the spectrum of the 
metric fluctuations and of the Bardeen potential. In fact, from Eq. f l4.14p the expression for 
B{x, t) becomes 



Mp a2(r 

Using Eq. f l4.24p and recalling that \E'(x, r 
is: 

ncBiw) 



a{T') q^{x,T') dr', Cb{w) 



3w(w + 11 



(4.24) 



-'HB{x,t) the mode expansion for \l/(x, r) 



^(f,r) 



Mp a2(r) (27r)3/2 
\{t')Ut') dr'. 



d^k 



(4.25) 
(4.26) 



The two-point function and the related power spectrum are simply 

sin A;r 



(^'(x,r)^'(y,r)) 

V^{k,T) 



d\nkVq,{k, t) 



kr 



\gk{T)f{2nk + l). 



(4.27) 



Within the same logic it is straightforward to derive the power spectrum of the density 
contrast. From Eqs. f l4.15p and fl4.16p after some algebra it can be shown that 



6t(x, t) 





0_ 




3Hcl\ ) 


dr 





(4.28) 
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Finally, using Eq. f l4.28p the power spectrum of the density contrast turns out to be 



hk(r) 



C|H|/ifc(r)|'(2rI, + l), 



2 7r2 M\ 



d 



H dr 



fkir) 



Cx(w) 



l3(w 



w 



a 



l+3c2 



(4.29) 



Let us now suppose that the initial fluid phase is dominated by thermal phonons. In this 
rather realistic situation w = 1/3 and Cs = y/w. From Eqs. (14.22^ and (14.23^ the spectrum 
of curvature perturbations can be recast in the following form: 



N, 



T 



1440^3 VA;/f/ VM 



eff (J^ 



coth 



eff 



360 



_k_ 



T 



(4.30) 



aH = "H; in Eq. f l4.30p the physical temperature T is related to the Hubble rate 
as T^/30 where Nf;ff denotes the effective number of relativistic degrees of 

freedom. In the limit {k/kH){T/Mp) < 1, Eq. (14.301) becomes: 



where kn 



Vnik, kn, T) 



eff 



k 



24 VSOtt Vfc/f / \M 



T 



(4.31) 



When the wavenumber is of the order of the particle horizon during the protoinflationary 
phase the amplitude of the curvature phonons is solely controlled by the temperature which 
must not exceed the Planck temperature. When k k^ the power spectrum is further 
suppressed. Since the transition to the fully developed inflationary phase occurs for r ~ 
r* ~ l/(a*iJ*) the spectrum computed from Eqs. (14.301) and (14.311) stops being valid for 
k ~ k^: = This means that {k/k^) < 1 implying k/ku < \'t/'t*\ which simply tells 

that the initial conditions during the protoinflationary phase must be set not too early or, 
equivalently, not too close to the Planck curvature scale. 

Similar considerations apply in the discussion the spectrum of the metric perturbations 
and of the density contrast. From Eqs. (I4.24l) - (l4.27p we arrive at the following explicit 
expression: 



V^{k, kn, T) 



V3N, 



eff 



360 



T 



coth 



eff 



360 



k 

kn 



(4.32) 



Unlike Vn^k, kn, T), Vif,{k, kn, T) is always a sharply decreasing function for k > kn- From 
Eq. (I4.29P the spectrum of the density contrast is: 



2 /k^ + 12kl 



H 



A;2 



coth 



eff 



360 



k 

kn 



T 



(4.33) 



where fc* = a*if* and, as before, kn = aH. From Eq. (I4.33P it can be argued that as long as 
k < k^: and H^, < Mp the modes inside the Hubble radius (i.e. k ^ kn) do not jeopardize 
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the validity of the perturbative expansion. During the protoinflationary phase the Hubble 
rate sharply increases towards the singularity and, in this situation, it can happen that large 
fluctuations arise for typical scales larger than the Hubble radius. This effect is caused by 
the nearness of the singularity and he lower limit in the time coordinate should be fixed by 
enforcing the validity of the perturbative expansion. 



4.3 Suppression of density contrast and metric fluctuations 



During the slow-roll phase and and in the limit Zpr — )■ the analog of Eqs. fl4.14p - fl4.16p can 
be written as 



B + 2HB- — 
a 

7^ = — , 

a e 

2V^B 
* ~ ^ Ha 



0, 



2e7^, 



St 



Sp<p 



(4.34) 
(4.35) 
(4.36) 



where 6^ is now dominated by the fluctuations of the inflaton. We shall assume, for instance, 
the validity of the solution f l4.4p with /3 ^ 1. This requirement is even too restrictive since 
the results discussed hereunder are valid in the slow-roll approximation, i.e. when both 
slow-roll parameter^ 



H 



Ml 



V \ 2 



V, 



(4.37) 



are much smaller than 1 but not necessarily constant. If TZ is continuous across the transition 
Eq. fl4.34p implies that the expression for B{x,t) becomes: 



7^(f,^') a(t') 



dt'. 



(4.38) 



Conversely, the continuity of B{x,t) in Eq. fl4.35p implies 7?. ~ for typical wavelengths 
larger than the Hubble radius. Therefore Eqs. fl4.16p . fl4.36p and fl4.38p imply 



B{X^ tpy.) 



B{X, tinf) 



TZ^{x) 
cipj- Hpj- 
7^*(x) 

Ctpr Hpj- 



a J 

1 / a(t') dt' 

a J 



7^*(x) 

O^inf Hinf V e + 1 



(4.39) 
(4.40) 



where tpr and t^f are, respectively, the values of the cosmic time coordinate during the 
protoinflationary phase and during the slow- roll phase when e ^ 1 and 1] <^ 1. Recalling 



^The slow-roll parameter e must not be confused with the parameter of the gauge transformation intro- 
duced in Eq. ()2.14|) . These two variables never appear together either in the preceding or in the following 
discussion so that no confusion is possible. 
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that = -{aH)B{x,t), Eqs. fO^]) and f OD]) show that 



^(^, w) = (ttt) ^^-^^^ 

Since during the protoinflationary phase the evolution is by definition decelerated, < w < 1. 
Conversely in the infiationary regime e ^ 1 which demonstrates the suppression of the 
metric fiuctuation. The same logic leads to following relation valid in the case of the density 
contrasts: 

St{x, tirrf) = ^^^"^^^^ 6t{x, tpr). (4.42) 

Once more, since 3{w + l)/2 ~ 0{1) and e <^ 1 we also have that 6t{x,tinf) ^ 6t{x,tpr) 
demonstrating the suppression of the density contrast across the protoinfiationary boundary. 
Note that the results reported here hold for a slow-roll parameter e which is not necessarily 
constant. 



5 Enhanced protoinflationary inhomogeneities 
5.1 Monotonicity properties 

The dynamics of the protoinfiationary transition depends on the behaviour of the pump 
fields Zpr{t) and z^(t). In the sudden treatment of the transition discussed in section H] the 
pump fields Zpr{t) and z^{t) evolve monotonically in cosmic time and the corresponding 
rates are positive definite (i.e. Zpr/zpr > and z^/z^, > 0). Furthermore the evolution of 
Zpr{t) is decelerated (i.e. Zpr(t) < 0). These requirements are satisfied in section H] where 
Zpr(t) oc a(t), the evolution is decelerated as long as the perfect barotropic fiuid dominates 
(i.e. Zprit) > and Zpr{t) < 0); furthermore, recalling Eq. (14.371) . during the slow-roll phase 



(p H 



H{l + ri + e) >0, (5.1) 



since, by definition of slow-roll, e ^ 1 and ?7 ^ 1 with e > 0. In the particular case of 
exponential potentials the slow-roll parameters are constant. The monotonicity requirements 
define the conventional protoinfiationary dynamics and they seem to be naively sufficient to 
implement a successful transition. Are they at all necessary? Is it possible to implement the 
transition in a different way without relying on the monotonicity of the pump fields? If this 
is the case, which are the consequences for the evolution of curvature perturbations? 

If the transition to the slow- roll phase is not sudden the behaviour of Zpr{t) and z^(t) is 
not necessarily monotonic. Whenever Zpr{t) and z^{t) do not evolve monotonically they can 
develop either a global maximum or a global minimum. There can be even more cumbersome 
evolutions where a finite number of maxima and minima arise. These situations can be 
discussed after having addressed the basic case where either Zpr(t) or z^{t) (or both) vanish 
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for a finite value of the cosmic time coordinate t. The absence of monotonic behaviour is 
dynamically reahzed in different ways by appropriately modifying the relative weight of the 
inflationary and fluid components in the total energy-momentum tensor of the system. A 
class of analytic solutions exhibiting non-monotonic behaviour for the pump fields can be 
obtained by solving Eqs. (I2.6p - (l2.9p whose explicit form, in the cosmic time coordinate, is: 



3H^Ml = ^ + Viip)+ppr, 



2 if Mp = -(f^ - {ppr + ppr) 



V3 + 3H(p + 



dV 



dip 

Ppr + 3H{ppr + Ppr 



0, 



0. 



(5.2) 
(5.3) 
(5.4) 
(5.5) 



We shall be interested in solutions where Zpr{t) and Zi^{t) are not monotonic but satisfy the 
correct boundary conditions typical of a protoinflationary dynamics. Equation (15. 3p can be 
written in terms of z^{t) and Zpr{t) and using the definition of e{t) given in Eq. (I4.37P 



a'MUit) = zi{t) + cyp^{t) 



(5.6) 



5.2 Exact solutions 

By requiring that Zpr{t) and z^p{t) are proportional it is possible to obtain a suitable ansatz 
for the solution of the whole system subjected to the requirement that the evolution of z^{t) 
and Zpr{t) is not monotonic. In the case of perfect barotropic fluid with constant sound 
speed the full solution of Eqs. fl2.6p - fl2.7p and (l2.8p -f l2T9|) can be expressed in terms of the 
scale factor and of the inflaton field: 



a{t) = a* 



sinh {(3 H, t) 



/3 



3(w + l) 



tanh 



2 



V 2 



(5.7) 
(5.8) 



where the parameter f2* and the protoinflationary energy density are defined, respectively, 
as 

p* 



Finally the inflaton potential is 



3Hm 



=2-' 



3H'^mI + -{1-w)H',mI{1 



P*\ — 

a 



fi*) sinh 



(5.9) 



1(3 ((^ - v^o) 



2 (l-fi,)Mp 



(5.10) 



Equations fl5.7p - fl5.10p solve Eqs. fl5.2p - fl5.5p in the case of a constant barotropic index. 
Furthermore, as anticipated, the solution satisfies the boundary conditions characterizing 
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the protoinflationary transition. In particular for l3H^,t < 1 the solution is decelerated and 
from Eq. fl5.7p we have 

ait)^a^i(3Hjy/^, H, = l-= ^ (5.11) 

pti 6[w + l)ti 

In the opposite limit (i.e. ^ 1) the solution is accelerated with H{t) ~ H^,: 

^^^^ = tanh(/3i7,t)' ^ ^ ~ smS m^t)' ^^'^^^ 

The parameter Vt^ = p*/(3if^Mp) < 1 measures the amount of radiation at the onset of 
inflation. From Eq. fl5.1ip . H^, is not the initial curvature scale but rather the curvature 
scale at the onset of the inflationary evolution. The conditions a > and a > imply, from 
Eq. dEB, 

a = ^^^[sinh (i7,t/3)]^/^"^| 1-2/3 + cosh [2i/,/5t] I > 0. (5.13) 

From Eq. flHTT^ a > iff cosh\/3 H,t) > (;9 - 1). The be ginning of the inflationary phase 
is then given by 

HJ^ = 1 In [^^ + y^^] . (5.14) 

This requirement shows that ranges between 0.54 if^^ in the case w = 1/3 and and 
0.49 in the case w = 1. For numerical purposes related to the evolution of the inhomo- 
geneities (see the discussion hereunder) it is practical to use H^, t as new time coordinate. 
In fact approximately coincides with the natural logarithm of the total number of in- 
flationary efolds. Hence -ff*tmax — C^(lniVtot) where we can take, for illustrative purposes 
A^tot between 70 and 100. Equations (14.371) and fl5.7p -( l578l) imply that e(t) gets progressively 
much smaller than 1 for P H^t ^ 1 while Tj{t) and //(t) are nearly constant in the same limit: 

/3 /X ^ /X ^ [1 + cosh^ (/3 i/, t)l 

e t = -, riit) = -13, nit) = ^ , * • 5.15 

cosh' (/3 i/, t) ' ' ' ^ ^ cosh" {(3 H,t) ^ ^ 

In conclusion the solution of Eqs. f l5.7p -( l5TT0l) satisfies the physical properties of a protoin- 
flationary regime and can be used to describe the transition to from a decelerated epoch to 
an accelerated phase where the slow roll condition is verified for e but not for rj. Notice, 
finally, that when fi,,, ^ 1 the potential V{(p) of Eq. fl5.10p can always be expanded in 
powers of i7* as 

viif) = Voiif) + v,iif)n, + V2iip)nl + ... (5.16) 

where the functions Vo{(f), Vi{(f), V2{f) are uniquely fixed from the exact form of the po- 
tential. 

The class of solutions reported in Eqs. (I5.7p - (l5.10p can be used to determine the explicit 
forms of Zpr{t) and of z^{t): 



/ 3(i + w) a{t) -jj Zprjt) _ 1 yn; 
'''■^^^ - V^i^^ *cosh(/3i7.t) ^-"T^Tr^- ^ ^ 
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Equations f l5.17p are illustrated in Fig. [T] where the non-monotonic behaviour of the pump 
fields is evident. In Fig. [T]the scale is linear on both axes and the three curves correspond to 
different choices of the barotropic indices. Since the plots of Fig. [T]are presented in terms of 
they hold for any value of H^. In spite of this the value of if* must be assigned in the 
numerical integration of the evolution of the inhomogeneities (see the discussion reported 
hereunder) since t controls, ultimately, the absolute normalization of the power spectra. As 
Eq. fl5.17p shows, the value of fi^, determines the relative magnitude of Zp^ and z^p. The 

Q. = 0.001 =0.001 




Figure 1: The evolution of Zpr{t) and of z^p{t) is illustrated in the case of the exact solution 
reported in Eqs. (15.71) and (15. 8p . 

difference in the overall sign of Zpr{t) and z^p{t) comes from the sign of (p which has been 
taken to be negative in Eq. (15. 8p . A flip in the sign of entails a flip in the sign of z^{t). 



5.3 Numerical integrations 

The system of Eqs. (I3.12p and (I3.13P will now be integrated in Fourier space and in the 
cosmic time parametrization which is preferable since the class of exact solutions reported 
in Eqs. (I5.7p - (l5.10p has a simpler expression in cosmic time. Equations (I3.12p and (I3.13P 
can be written in the form of a plane autonomous system: 

Qv = Pv, (5.18) 

Qx = Px^ (5-19) 

Pv = -Hp^ + A^^{k,t)qy + B^^(t)q^ + Cy^(t)p^, (5.20) 

Px = -Hp^ + A^^{k,t)q^ + B^^{t)q^ + C-^^{t)p^. (5.21) 

Equations (I5.18p -( l5.2ip hold in Fourier space, i.e. g„ = qv{k,t), q^ = qy^ik^t) and similarly 
for p„ and p^^. For notational convenience and to avoid the proliferation of subscripts the 
reference to the modulus of the wavenumber k has been omitted unless strictly necessary. 
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The coefficients A^y{k,t), By^{t) and Cy^{t) are given, respectively, by: 

Ayy(^k,t) - 



Sw-1- (3W- l)(3w + l) c^P 
H + -H^ - 



+ 



By^it) 



2 

P + P 

\/Ppr ~\~ Ppr i ( V 



4 



w - l)^"^ - 2 3(w + l)H^ + 2H 



dV 



w - 1)0^ - 2 3(w + 1)H^ + 2H 



- 2{w + l)— + 2{l-w)Hip 
dip 



The coefficients A^^{k,t), B^y{t) and C^i>(t) are instead given by: 



AH 



+ [jj) + ^H') + {p,r + Ppr) 

y/Ppr + Ppr \ (w - l)(3w - 1) ^ £V_ 

w dip 



4/7 



4(if + m^) + (pp, + Pp.) 



w — 1 



w 



w — 1 



w 



w 



^Vpr + Pp r fW - 1 

2H 



cp. 



(5.22) 

(5.23) 
(5.24) 



(5.25) 

(5.26) 
(5.27) 



Equations (ESS])- (ESI) and Eqs. (15:251) - ( 15:271) are related to Eqs. ( 137[4D -( 13J6D and (13X71) - 




Figure 2: The spectrum of g^^ (plot at the left) and of (plot at the right). Common 
logarithms are of the corresponding quantities are used on both axes. 
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fl3.19p modulo an overall sign difference! and the obvious algebraic changes related to the 
differences between the conformal and the cosmic time parametrization. Notice, furthermore, 
that the coefficients A^^^k, t) and A^^{k, r) depend both on the cosmic time coordinate and 
on the wavenumber which comes from the three-dimensional Laplacian of and g-^. 

The initial conditions for the numerical integration of Eqs. f l5.18p -( 15.2ip are fixed by 
requiring that, at the initial integration time, Qx ~ Px ~ ^ while qy and are determined 
from the thermal phonons during a radiation-dominated regime (i.e. w = 1/3). The 
considerations of section [Hand the exact solutions presented in this section allow for different 
values of the barotropic indices and for a wider set of initial conditions. In spite of this it 
seems less confusing, for the purposes of the presentation, not to indulge in an excessive 
multiplication of examples which could cast shadow over the main motivation of the present 
analysis. 

Since the initial conditions are given during the protoinflationary phase (i.e. H^,t <^ 1), 
nearly all the relevant modes are larger than the particle horizon since k/a{ti) ^ H{ti). This 
aspect can be understood by looking at the sharp increase of H{t) in the limit t — i- where 
the Hubble rate dominates over l/a{t) in spite of the possible largeness of the wavenumber. 
The full, dashed and dot-dashed curves appearing in the plots of Figs. [2] and [3] correspond 
to different values of the wavenumbers k which are expressed in units of in terms of the 
rescaled wavenumber k = k/H^ (see also the legends of each plot). The values of if^, used 
in each numerical integration are reported, in natural Planckian units Mp = 1, in the title 
of each plot illustrated in Figs. |2] and [31 

As explained in sections [3] and HI once the spectra of and q^ are known, all the other 
quantities can be computed such as the spectra of B and or the spectra of \E' and TZ. On 
both axes of Figs. [2] and [3] the common logarithm of the corresponding quantity is reported. 
Instead of Uustrating the power spectrum it is more practical to compute the square root of 
it, as indicated on the vertical axes of Figs. |2]and|3l During the protoinflationary phase (i.e. 
roughly H^t < 1) the spectrum goes as k^/^ as it can be deduced from the different curves 
holding for decreasing values of k. Recall that in Figs. [2] and [3] the common logarithm of 
the corresponding quantity is reported on both axes. The slope k^/^ arises since, initially, 
k/{2Ta) <^ 1 and therefore V^J^ oc k^Ta/k. As H^t ^ 1 the spectrum tends to change 
even if this aspect is more visible by looking directly at the spectrum of TZ reported in Fig. 
|3l From Fig. |3]the spectrum of curvature perturbations exhibits a sudden growth which can 
even be understood on a qualitative basis. Recall, in this respect, the general expression of 
TZit) in terms of q^ and q^ (see Eq. f l3.2ip ). When the spectrum of q^ and q^ is bounded (see 
Fig. [2]) Vn increases when either Zpr{t) or Zip{t) decrease. The slope of the spectrum goes as 
fQj^ ^ \ while it is quasi-flat for H^t > 1. In the left plot of Fig. [3] the parameters 

and if* coincide with those adopted in Fig. |2l in the right plot of Fig. [3] the values of 
the fiducial parameters are instead fi* = 10~^ and if* = 10~^. The examples discussed here 

^Note that in Eqs. (|3.12p and (|3.13p the coefficients all appear at the left hand side in the corresponding 
equations. In Eqs. (|5.20p and (|5.2ip the coefficients appear at the right hand side. 
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w=1 /3, il. =0.001 , H. =1 ^ w=1 /3, n. =0.00001 , H. =1 0"' 




Figure 3: The spectra of the curvature perturbations on comoving orthogonal hypersurfaces 
for shghtly different sets of fiducial parameters. As in Fig. |2] common logarithms of the 
corresponding quantity are used on both axes. 

are just meant to show that the lack of monotonicity of Zpr{t) and Zip{t) always entails a 
potential decrease of the pump fields and, hence, a potential increase of the total curvature 
perturbations. It is interesting to point out that the possibility of anomalous growth of 
curvature perturbations in single field inflationary models has been discussed some time 
ago but within a different perspective [27] and not in connection with the protoinfiationary 
dynamics. 

The analysis of this section has been conducted in the case of a single protoinfiationary 
fluid and a single inflaton; there seems to be however no obstruction for the possible general- 
ization of this approach to a finite number of perfect fluids simultaneously present together 
with a finite number of inflaton fields. Finally, the considerations of this paper just assumed 
the validity of the perturbative expansion prior to the onset of the inflationary evolution. 
The overall consistency of the approach, as already stressed in section HI implies that valid 
conclusions can only be drawn in the regions of the parameter space where the perturbative 
expansion is valid. At the same time there are techniques, like the gradient expansion, which 
can be apphed in situations where the conventional perturbative expansion fails [281 129| 130]. 
This analysis is anyway beyond the scopes of the present paper. 

6 Concluding remarks 

The effects of a dynamical phase preceding inflation are customarily parametrized in terms 
of an initial number of inflaton quanta modifying the standard vacuum initial conditions for 
the evolution of the large-scale inhomogeneities. In this paper the implications of a protoin- 
fiationary phase of decelerated expansion have been examined with the purpose of relaxing 
the standard lore. The present approach stipulates that, initially, the energy-momentum 
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tensor is dominated by an irrotational fluid. After a transient regime (wliicli may be eitlier 
sliarp or delayed) the inflaton potential dominates the total energy-momentum tensor and 
the slow-roll dynamics starts off. The initial conditions for the large-scale inhomogeneities 
during the protoinfiationary phase can be either quantum mechanical or, more realistically, 
thermal if the protoinfiationary phase is dominated by radiation. 

Is the existence of a protoinfiationary phase just equivalent to a modified set of initial 
conditions manually imposed on the modes of the perturbations? The results of the present 
analysis show that the answer this question is positive insofar as the protoinfiationary tran- 
sition is sufficiently sharp and satisfies a number of monotonicity requirements. Across the 
protoinfiationary boundary it is plausible to demand that the extrinsic curvature is contin- 
uous while the stress tensor undergoes a finite discontinuity on the constant energy density 
hypersurface. In this picture the transition is sharp and both the density contrast and 
the metric fluctuations are suppressed across the transition. This sudden approximation 
merely extends to the protoinflationary boundary the standard discussion of the various 
post-inflationary transitions. 

The dynamics of the inhomogeneities is determined by the evolution of the two quasi- 
normal modes of the system. Depending on the time evolution of the pump fields controlling 
the contribution of the quasi-normal modes to the curvature perturbations, the protoinfia- 
tionary transition is naturally classified in two broad classes. In the first category the pump 
fields are monotonic and the quasi-normal modes of the system are nearly not interacting. 
The second category contemplates the complementary situation where the transition is not 
monotonic and the pump fields have, at least, either a maximum or a minimum. In the lat- 
ter case large-scale curvature perturbations can be enhanced in comparison with the results 
obtainable in the case of monotonic transitions. 

The large-scale modifications of the temperature and polarization power spectra are often 
engineered by assuming that the onset of infiationary dynamics is only characterized by a 
single infiaton field possibly with non-standard initial conditions for its inhomogeneities. 
The exact solutions obtained in this paper, as well and the numerical treatment of the 
corresponding large-scale inhomogeneities, suggest that the protoinfiationary dynamics is 
richer than expected. In a more conservative perspective the obtained results show that the 
potential enhancement of large-scale curvature perturbations can be avoided by demanding 
that the evolution of the pump fields is strictly monotonic. It is however unclear whether 
or not the monotonicity of the pump flelds must be considered a generic feature of the 
protoinflationary dynamics. 
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A Explicit derivation of the quasi-normal modes 



In this appendix the details of the derivations leading to Eqs. fl3.12p and fl3.13p will be 
outlined. By taking the first time derivative of Eq. f l2.33p the following equation arises: 
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u" + H'il - ?>w)u + nil - ?>w)u 



w + l 



5' + 



(A.l) 



where, as already mentioned in section [3l the barotropic index w has been taken to be 
constant so that Cg = a/w; note that 5 = Sppr/ppr denotes, only in this appendix, the density 
contrast of the protoinfiationary fluid. Now the logic is to express the right hand side of Eq. 
f lA.ip solely in terms of u and x- The same procedure must then be applied in the case of 
Eq. i^M>- 

The right hand side of Eq. f lA.ip is expressible in terms of u, x and x'- To achieve this 
goal, Eqs. (I2.22P and f l2.32p can be summed up term by term. For a perfect barotropic fluid 
Eq. f l2.32p is given by 5' = — (1 + w)Opr] by then using the expression for 0' obtainable from 
Eq. f l2.29p it is easy to show, after some algebra, that 
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Equation (13. lip will then be used into Eq. (]A.2p to eliminate and, therefore, the final form 
of Eq. CO) is: 
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The same strategy used to derive Eq. (]A.3P leads to the evolution equation of the second 
quasinormal mode associated with Eq. (12.3ip . Here the problem is to trade the last three 
terms appearing at the right hand side of Eq. (12.3 ip for appropriate combinations of u 
and u'. The relevant combination arising from the terms containing 0, V^B and 
transformed by using Eqs. (I2.22p . (I2.29p and (13.110 . The final result is 
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The derivation of Eqs. flA.Sp and f lA.4p breaks down when Cg = y/w = so that the latter 
case must be separately treated. The starting point of the derivation will be the same except 
that, since 5ppr = 0, Eq. f l2.33p is even simpler. Repeating all the steps of the derivation we 
obtain, after some algebra, that the evolution equation for u becomes: 
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Similarly, in the case Cg 



w = 0, the evolution equation for x becomes, after some algebra: 
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Going back to the main derivation and to Eqs. (1A.3P and (]A.4p . the expression for the 
evolution of u can be modified by defining a new variable v such that 
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In Eq. (1A.7P Planck units Mp = 1 have been set and, in these units, the equations for v and 
X become, respectively. 
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From Eqs. flA.Sp and f lA.Qp the evolution equations for the canonical normal modes 

av 

qv = — , gx = «x, (A.io) 

leads to Eqs. fl3.12p and f l3.13p . The variable v is the closest analog, in the fluid case, to 
the infiaton fluctuation x- If is useful to recall, as a technical remark, that the intermediate 
steps involving the variable v can be avoided. In the latter case is directly expressible in 
terms of the variable T-L u; with this strategy the algebraic derivation is, though, 

more tedious. 
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